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ABSTRACT 



With the possibilities of high spatial resolution imaging and spectroscopy as well as infrared (IR) interferometry, the dusty environ- 
ments (= "dusty torus") of active galactic nuclei (AGN) are now in reach of observations. Following our paper I on ground-based 
mid-IR spectro-photometry (Honig et al. 2010), we present an upgrade to our radiative transfer model of 3-dimensional clumpy dust 
tori. The upgrade with respect to Honig et al. (2006) concerns an improved handling of the diffuse radiation field in the torus which 
is approximated by a statistical approach. The models are presented as tools to translate classical and interferometric observations 
into characteristic properties of the dust distribution. We compare model SEDs for different chemical and grain-size compositions 
of the dust and find that clouds with standard ISM dust and optical depth Tv ~ 50 appear in overall agreement with observed IR 
SEDs. By studying parameter dependencies, it is shown that type 1 AGN SEDs, in particular the mid-IR spectral index, can be used 
to constrain the radial dust cloud distribution power-law index a, while other parameters are more difficult to assess using SEDs only. 
Interferometry adds important additional information for modeling when interpreted simultaneously with the SED. Although type 2 
AGN can, in principle, be used to constrain model parameters as well, obscuration effects make the analysis more ambiguous. We 
propose a simple, interferometry-based method to distinguish between "compact" and "extended" radial dust distributions without 
detailed modeling of the data and introduce a way to easily determine individual or sample average model parameters using the 
observed optical depth in the silicate feature and the mid-IR spectral index. 
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1. Introduction 

The Unification Scheme of active galactic nuclei (AGN) is re- 
lying on a toroidal region filled with molecular gas and dust to 
, explain the observed dichotomy of broad- and narrow-line AGN 

■ (lAntonuccilll99llUrrv & Padovanilll995h . The presence of this 
" optically-thick region is well supported by direct and indirect 

observational evidence. Miller & Antonucci ( 1983) showed that 
the narrow emission-line galaxy NGC 1068 shows broad op- 
, tical emission lines in polarized light, scattered outward from 

■ the innermost region which is obscured by the "dust torus". 
' Several succeeding studies confirmed the presence of a hidden 
, broa d-line AGN in many narrow-line objects (e.g. Moran et al. 
' 12000 , and references therein). From number statistics of ob- 
' scured type-2 AGN and non-obscured type-1 AGN as well 

as from direct observations of the pc-scaled molecular gas in 
nearby AGN, it is inferred that the torus is not only optically 
thick but also geometrically thick, with a cha racteristic scale 
height h/r ~ 1 (e.g. [M aiolino & Rieke 1995; Lacv et al.ll2004t 
lMartuiez-Sansigre etal.ii2006i;. Hicks et al...2009'!K 

Obscuration of the torus does not only affect the optical 
broad lines but can also be observed in X-ray and optical/UV 
continuum emission from the accretion disk and its immedi- 
ate vicinity. Optically-obscured type-2 AGN are usually asso- 
ciated with high Hydrogen column densities in the X-ray, from 
~ 10^^ to Compton-thick columns of > 10^^"^^ cm"^ (e.g. 
IShietal.ll2006b . The material causing these high column den- 
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sities is generally believed to reside in the dust torus, tracing 
the gaseous component which dominates the total mass of the 
torus. Moreover, strong reflection on parsec scales, as observed 
in the X-ray continuum around 30keV and associated with large 
equivalent width of the Fe-Kor line, are interpreted as signs for 
the presence of the torus. 

Since a sizable portion of the torus consists of dust, the 
UV/optical accretion disk emission heats up the dust which ther- 
mally re-radiates the received energy in the infrared at temper- 
atures of some lOOK up to the dust sublimation temperature of 
~1500K. These temperatures are reached on scales of about 0.1 
to lOOpc, depending on the actual AGN luminosity, torus ge- 
ometry, and dust composition. The thermal radiation from the 
torus dominates the total near- to mid-infrared (IR) emission of 
the AGN ("red bump") with a noticeable cut-off at around 1 //m 
where the UV/optical accretion disk emission starts to dominate 
("big blue bump"). At around 10 and 18 //m, type-1 AGN show 
some broad spectral features caused by hot silicate dust which is 
believed to be associated with the inner part of the torus. On the 
other hand, these silicate features appear in absorption in most 
type 2 AGN where only cooler dust is seen. 

Due to the small spatial scales at which the torus resides, 
it is a difficult task to directly resolve it by infrared observa- 
tions. Thus, first reports of successful resolutio n of the nucleus of 
an AGN involved interferometric techniques. IWittkowski et al.l 
(119980 used bispectrum speckle interferometry to resolve the 
nucleus of the Seyfert 2 galaxy NGC 1068 in the /T-band, 
foUowed-up by additiona l //-band observations dWeigelt et al.l 
l2004 . ISwain et al. ] (l2003h presented /T-band long-baseline inter- 
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ferometry of the type-1 AGN NGC 4151 which arguably resolve 
the in nermost hot region of the dust torus (see Kishimot o et al.l 
120071 for this interpretation). More recently, VLTI/MIDI long- 
baseline mid-IR spectro-interferometry directly revealed the 
parsec-scaled dust emission sizes in the wavelength band from 
8 to 13^m for a riumber of nearby type 1 and type 2 



ivt: 

AGN (Jaffe et al. 2004: 'Tristram et al."2007:'Becker t et alj2008t 
[Raban et al. 2009; Tristram et al. 2009). These observations fi- 
nally confirmed the basic picture of the dust torus while more 
detailed characteristics remain unclear. 

The dust torus has been the subject of several kinds of 
models in order to extract physical properties from broad SED 
and, in some cases, interferometric observations. Initially, most 
authors used smooth dust distributions with different kind of 



radial and vertical dens i ty profiles (e.g. Pier & Krolik 1993 
Gran ato & Daneselll994 l Efstathiou & Rowan -Robinson 1995 . 
Sch artmann et al. 2005; Frit z et al .II2006^ . t was, however, early 
noted that the dust is most probably arranged in clouds in- 
stead of bein g smoothly distributed (e.g. Krolik & Begelman 
ri988l; Itacconi et all Il994i; iHonig & BeckertI 12007). This idea 
received further support by interferometric observations which 
arguably rule out smooth dust distributions dJafFe et al.l l2004t 
iTristram et al] |2007). Several radiative transfer models have 
been developed to account for 2D or 3D clumpy dust distri- 
butions (' Nenkova et al. 2002; Dullemond & van Bemmelll200l 
iHoniget al...2006; .Schartmann et al...2008l) . All of these mod- 
els appear in more or less good agreement with observations, 
while the res ulting torus or cloud properties differ signifi- 
cantly. iHonig e t al. (2006) used optically thick dust clouds and 
a low torus volume filling factor to simultaneously model near- 
and mid-IR photometry and i nterferomet ry of NGC 10 6 8 (see 
also IHonig et"al] l2007i lIOOSi) . while Sch artmann etaTI ( l2008l) 
model a torus with high volume filling factor and optically thin 
clouds for a similar set of observations of the Circinus galaxy. 
Although the actual properties of the dust clouds are not yet con- 
strained by observations, theoretical predictions and hydrody- 
namic simulation s arguably favor small and compact optically- 
thick clouds (e.g IV ollmer et al."2004'; 'Beckert & Duschll 120041: 
iHonig & Beckertll2007i: .Schartmann et al. 2009ll . 

In this paper, we present an update of our radiative t ransfe r 
model of 3D clumpy AGN tori presented in lHonig et alj (l2006h . 
In particular, a better handling of the diffuse radiation field in- 
side the torus has been implemente d to overcome the limitations 
discussed in the Appendix A. 2 of 'Honig et al.l (1200 6). and the 
possibility of different dust compositions and grain sizes is now 
included. The impro ved handhng of the diffuse radiation field 
has been proposed in lHonig I (1200 8) and is described in more de- 
tail in Sect. 12.3! We intend to provide models for observers and 
come up with some simple correlations between observational 
and model parameters. Some emphasis will be put on the radial 
distribution of the dust and dust grain size and chemical com- 
position, the latter being only marginally explored in literature 
despite observati onal evidence that the grain composition might 
be crucial re.g. lKishimoto et al.ll2007ll2()09al) . 

In Sect.|2]we outline our model strategy and methods and in- 
troduce the mathematical background of the clumpy torus model 
including model parameters. Following in Sect. [3] simulation re- 
sults of individual dust clouds are presented. These simulations 
build up the dust cloud databases which serve as an input to the 
torus model simulations. In Sect. |4] we show our results from 
the torus simulations and discuss various aspects of interpret- 
ing SEDs and IR interferometry. The results are summarized in 
Sect.H 



2. Torus model fundamentals 

In this section, we will describe the mathematical and strate- 
gic base of our torus model and discuss underlying assumptions. 
First we outline the general model strategy and compare it con- 
ceptually to torus models in literature. Then, we describe the 
theoretical basis and parameters involved in the model. 



2.1. Model Strategy 

The most straight-forward way of modeling clumpy dust tori in 3 
dimensions is direct Monte Carlo simulation of a geometrically 
well-defined, statistically arbitrary distributi on of dust clouds 
spread on a model grid, as demonstrated by ISchartmann et al.l 
(I2OO8I) . There are, however, several practical problems that oc- 
cur. When solving the radiative transfer equation by Monte Carlo 
simulations, it is important to properly sample optically thick 
surface regions by enough grid cells so that each cell is opti- 
cally thin. Otherwise, emission temperatures will be underes- 
timated leading to a wrong final source function which affects 
torus SEDs and images. In principle, adaptive grids can be used, 
but this may become a difficult task when aiming for > lO"' - lO'* 
randomly-arranged clouds. Thus, the number of model clouds 
has to be small, the total optical depth of each cloud has to be 
limited, and/or the volume filling factors have to be rather large, 
of the order of Oy ~ 1 (see also Dullemond & van BemmeJ 
I2OO5I) . ISchartmann et all (l2008l) showed that the SED of the nu- 
cleus of the Circinus galaxy as well as the position-angle- and 
baseline-dependence of the visibility is in good agreement with 
this kind of Monte Carlo models. Still, such direct Monte Carlo 
simulations take a lot of time and are, thus, not very flexible for 
modehng of observations. 

As mentioned in Sect.[T] it is most likely that individual torus 
clouds are actually optically thick (see also Sect. 13.1b . Based 
on optically thick clouds and the assum ption of a low volume 
filling factor, (i>v <K 1, iNenkova et al.l (I2OOZ) used a proba- 
bilistic approach for cloud heating and obscuration, depending 
on several model parameters (see also 'Natta & Panagial 119841 ; 
Nenkova et al. 2008a). Nenkova et al. (2008b) presented model 
SEDs simulated via this approach whic h are in general agree- 
ment with observations. iBeckertI (|2005|) used a similar proba- 
bilistic model to reproduce the high spatial resolution SED of 
NGC 1068. The benefit of this modeling approach is its time- 
efficient calculation of average model SEDs and average bright- 
ness distributions. On the other hand, as noted bv lNenkova et al.l 
(2008b), the probabilistic method is not capable of predict- 
ing variations in the overall SED and brightness distribution or 
small-scale spatial and position-angle variations (="clumpiness 
variations") of interferometric visibilities and phases caused by 
the random distribution of dust clouds in the torus. 

In iHonig et'al] (l2006h we presented a method to time- 
efficiently simulate torus model images and SEDs and account 
for the 3D statisti cal nature of clumpy dust distributions. The 
strategy follows ,Nenkova et al. I ([2008 a), however resolving the 
need for a probabilistic treatment of the dust distribution. Just 
as in the probabilistic approach, our model strategy is based on 
separating the simulations of individual cloud SEDs and im- 
ages, and the final SEDs and images of the torus (see Fig. [1] 
for an illustration of our method). In a first step, the phase- 
angle-dependent emission for each cloud is simulated by Monte 
Carlo radiative transfer simulations. For that, we apply the non- 
iterative method delineated by Biorkman & Wood (2001), based 
on work by Lucy (1999). In our previous model, we used the 
code mcSim (^Ohnaka et al...2006 ) which is able to simulate a 
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Fig. 1. Flow chart of our method for radiative transfer simula- 
tions o f 3D clumpy tori. The strategy follows Nenkova et al. 
(I2008ah . however resolving the need for a probabilistic treatment 
of the dust distribution. In a first step, we simulate databases of 
directly- (see Sect. 13.1b and indirectly-illuminated clouds (see 
Sect. 13.3.2b . The step using the indirectly-illuminated clouds as 
a heating source for directly-illuminated clouds (dashed arrow) 
would close an iterative loop which we argue in Sect. |23] con- 
verges satisfactorily after the initial loop. In the second step, 
clouds are randomly distributed in a torus according to the model 
parameters (see Sect. 12.2b . According to their local environ- 
ments, the distributed clouds are associated with database entries 
(see Sect. 12.3b . The emission from each cloud is then followed 
through the torus and final torus SEDs and images are obtained. 



variety of geometries. However, in order to obtain fast results 
for different dust compositions, we created a new Monte Carlo 
code which is optimized for the AGN-cloud-configuration. For 
each given dust composition, we determine the sublimation ra- 
dius Tsub - '"(7'sub - 1 500 K) from the source, and simulate sev- 
eral clouds at different distances (normalized for rsub)- To finally 
simulate the torus emission, dust clouds are randomly distributed 
around an AGN according to some physical and geometrical pa- 
rameters. Each cloud is associated with a pre-simulated model 
cloud, after accounting for the individual cloud's direct and in- 
direct heating balance (see Sect. 12.31 ). The final torus image and 
SED is calculated via raytracing along the line-of-sight from 
each cloud to the observer. Contrary to the probabilistic clumpy 
torus models, this method accounts for the actual 3-dimensional 
distribution of clouds in the torus involving all statistical varia- 
tions of randomly distributed clouds in a time-efficient way. 



2.2. Torus parameters 

The distribution of clouds in the torus is initially characterized 
by six model parameters. The torus parameters are (1) the radial 
dust-cloud distribution power-law index a, (2) the scale height 
h or the half-opening angle Oq, (3) the number of clouds along 
an equatorial line-of-sight A^o, (4) the cloud radius at the subli- 
mation radius Ra-o in units of rsub, (5) the cloud size distribution 
power-law index b, and (6) the outer torus radius Ro^t- From the 
model parameters, several other physical properties can be de- 
rived. However, we will show that some of these parameters can 
be considered only marginally influential on the torus SEDs and 
images, so that less parameters need to be constrained by obser- 
vations. In the following, we will briefly show how the model 
parameters are associated with the dust distribution of the torus. 

The most fundamental parameters deal with the geometri- 
cal distribution of the clouds in the torus. We separate the 3- 
dimensional dust distribution into a radial distribution power- 
law, rjr oc (r//-sub)" and a vertical distribution r]-/g. For the vertical 
distribution, it has become common practice to use a Gaussian 
distribution to reproduce the smooth transition from type 1 to 
type 2 viewing angles as implied by observations. However, two 
different ways of defining this Gaussian have evolved from re- 
cent torus models. One possibility is to distribute the clouds 
perpendicular to the equatorial plane, i.e. in z-direction in a 
cylindrical coordinate system. The resulting distribution func- 
tion T]- OC exp(-z^/2//^) depends on the scale height H{r) - h - r 
at radial distance r. This implies that the torus flares with con- 
stant H{r)/r, as suggested by isothermal disks. Here, h is the fun- 
damental parameter describing the vertical distribution. Another 
way of distributing clouds can be best understood in a spherical 
coordinate system: instead of defining the distribution function 
perpendicular to the equatorial plane, it is also possible to dis- 
tribute the clouds along an altitudinal path, i.e. in spherical 6- 
direction with rjg oc exp(-6 /0q ). Here, 0o - (tt/I - 0o) is the 
half-covering angle of the torus which is the complementary of 
the half-opening angle Oq (see remarks in Table [TJ. Most of the 
simulations shown in the later sections will be based on the 6- 
distribution, but we will include results using the z-distribution 
when drawing final conclusions from the SED models. Table [1] 
shows the full distribution functions J]r, and ij^/g. The complete 
dust cloud distribution function 

T](r,z/e)^27TT],.-T]-J9 (1) 

is normalized so that Jj^°"' rjrdr - 1 and rj^/g = 1 for z = or 
= 0. By this definition, rir can be understood as the normalized 
number of clouds per unit length, so that 



NoJ]{r,z/0)=nRliPN 



(2) 



describes the (actual) number of clouds per unit length, where 
Pn denotes the cloud number density per unit volume. Note that 
r, z, H, and R^-fl are given in units of rsub- 

After having fixed the distribution of clouds, the total obscu- 
ration (or dust mass) of the torus has to be defined. The most 
convenient way to do this is by defining the mean number of 
clouds A^o that intersect the line-of-sight along a radial path in the 
equatorial plane. From A^o it is easy to calculate the total number 
of clouds that have to be randomly distributed in the torus. 



Mot = A^o 



r](r,z/e) 



dV 



(3) 



whi ch is essentially the integrated version of eq. (I2) (see 
also iHonig et"aD l2006t iNenkova et al.1 l2008ah . Here, R^i{r) = 
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Table 1. Cloud distribution functions and relations to physical 
parameters of the torus 



Physical property 



Function 



radial dust distribution 



'''' " 

vertical dust distribution 77- = exp (-z'/2H^^ 



dust cloud radius 
total number of clouds 



or;;,, = exp(-6'V^o^) 
NL 



Here, A^o and h are scaling constants which do not depend on 
the radial distribution. On the other hand, we can include a first- 
_order approximation of obscuration effects, which accounts for 
~a decreasing number of directly -illuminated clouds with radius, 
"by multiplying cr^ir) with exp i-N{r, z)) as described in the next 
Sect. 12.31 Finally, to obtain SEDs and brightness distributions in 
this simplified type-1 model, the luminosity lJ,'™' is calculated 
by multiplying the surface filling factor with the source function 
of the clouds S y{r) and integrating from rjub to Rout, 



3+0-2* 
n 2 , 



<^v= jNoJ]{r,z/d)Ra 



(4) 



If restricting to the equatorial plane, the expression (i>v;o(r) - 
<i>v{f',zoi9 - 0) is derived as shown in Table[T] 

We note that while consistency of our modeling approach 
requires <t>vir,z/0),< 1 V(r, z) or (r, 0), the torus brightness dis- 
tribution and SED is not explicitly depending on but rather 
on the combination of the individual parameters. Actually, ap- 
proximately the same overall SEDs and images can be obtained 
when reducing <^v by reducing Rdfi, leaving all other parame- 
ters fixed (and calculate A^tot clouds according to the new ^ci;o)- 
This can be illustrated by introducing a simplifi ed version of the 
clumpy torus model for type 1 AGN (see also iKishimoto et al.l 
[2009 a). In a face-on view of the torus, the clouds are projected 
onto a ring ranging from rsub to Ro^t- First, let us assume that 
most of the clouds are directly heated by the AGN. Then, a "sur- 
face filling factor" crs(r) can be described using cloud number 
density projected z-directiorQ If we consider 77,, the surface 
filling factor becomes 



nRl ■ PA, dz = {2nf'^ Noh 77, ■ r.^h 

00 



(5) 



with the explicit dependence crs(r) oc r". In this definition, the 
surface filling factor can be considered as a weighting factor of 
how much the clouds at radius r contribute to the total intensity. 



' Actually, (Ts(r] can grow to larger than 1 which corresponds to the 
case that projection caused clouds to overlap. However, crs(r) > 1 is 
only valid if the dust column is optically thin at that wavelength. 



r (tor) 



= 27r 



with/e„ = V2eoexp(-(?o'/4) Erf 2^ 
volume filling factor at z = Ov o = 4 NoRa-o -4t^ ■ r"*'' 

: i ' ^"'.1 

— Notes: ' r, z, H, and Rao are given in units of rsub; ^ Model pa- 
rameters: (1) a: radial distribution power law index, (2) h: scale height, 
h = H(r)/r = const, or half-opening angle 60, (3) Nq: mean number 
of clouds along equatorial line-of-sight, (4) R^ho'- cloud radius at r^ub (in 
units of Tsub), (5) b: cloud size distribution index, (6) Rout'- the outer torus 
radius; ^ We define the altitudinal angle = in the torus mid-plane. 
Therefore we use, for convenience, the half-covering angle 60 = 
instead of the half-opening angle So in the equations. 



^ci;0 '"sub ■ denotes the radius of a dust clouds at radial coor- 
dinate r (in units of rjub), with parameters /?ci;0 and b as defined 
above. In Table [T] the explicit expression for A^tot is given for 
both vertical z- and ^-distributions. 

Since we deal with a clumpy torus, it is most sensible to have 
a volume filling factor <i>v < 1 all over the torus. The volume fill- 
ing factor Oy can be calculated by multiplying the cloud number 
density per unit volume and the cloud volume Vd - 4/37rR^-^, 
so that 



r 



o-s(r) ■ Syir)rdr 



(6) 



Interestingly, L^'™' does not depend on R^i or b, so that it is dif- 
-ficult to constrain these parameters by observations. Instead, for 
the purpose of SED modeling, they should be selected in a way 
to fulfill the "clumpy criterion" < 1 . 

However, since R^i and b control the total number of clouds 
A^tot, they have influence on small scale surface brightness vari- 
ations resulting in slightly different SEDs for different random 
arrangements of clouds (see lHonig et alJ2006l) or position-angle 
variations, which might be interesting for interferometry (see 
Sect. |42] ). Aside from R^i-o and b, in Sect. 14.1.41 we wifl fur- 
ther show that /?out cannot be considered as a free parameter but 
has to be chosen in a sensible way depending on a. 

In summary, the key torus model parameters which are di- 
rectly accessible from observations are a, No, and h or Gq. In 
addition, the cloud dust properties (dust composition and optical 
depth) will have some influence on the overall torus emission. 

2.3. The diffuse radiation field in the torus 



As described in Sect. 12.11 our model strategy makes use of pre- 
calculated databases of dust clouds. The main database consists 
of clouds which are directly heated by the AGN. However, due 
to obscuration effects within the torus, some clouds will not be 
directly exposed to the AGN radiation since the line-of-sight to 
the AGN is blocked by other clouds. Such clouds are heated 
only by the emission from other clouds in their vicinity. These 
nearby clouds form a "diffuse radiation field" where mainly 
those clouds contribute to which are directly AGN-heated. In 
Honig et al. (2006), we used an upper-limit-approximation for 
the diffuse radiation which led to overestimation of the emission 
longward of 20 jum. Here we describe a way to recover the dif- 
fuse radiation field statistically. 

For given distance from the AGN and optical depth of each 
cloud, the parameters which determine the temperature and 
emission of a cloud are (1) the fractional cloud area, /dh, which 
is directly heated by the AGN, and (2) the fraction of directly- 
heated clouds in the cloud's vicinity, fiuir). /dh determines the 
strength of direct heating. For /dh - 1, the directly -heated 
re-emission of the cloud is equal to the emission of the pre- 
calculated database cloud at the same distance tdb = '"ci from 
the AGN. If /dh < 1, less direct energy is received and the 
cloud emission corresponds to a database cloud at distance tdb = 
fell -\//dh- The parameter fm(r) determines the energy that is 
contained in the diffuse radiation field. It can be approximated 
by the probability P()(r, z) that a given cloud at r and z (or ff) is di- 
rectly illuminated. Note that P^ir, z) can be considered a global, 
probabilistic version of /dh, which i s a local and individua l 
property of each cloud. According to iNatta & Panagial ( |1984|) . 
P()(r,z) can be derived from Poisson statistics as Po(r,z) = 
exp (-N(r, z/0)), where N(r, z/0) - Nq T](r, z/0)ds is the mean 
number of clouds along the path s from the center to (r,z/6). In 
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0.01 



0.001 




directly heated clouds ig„ 
not directly heated clouds (l-fpH) 
diffuse radiation field ^luO-^au) 



distance (r^^^) 

Fig. 2. Illustration of the contribution of directly- and indirectly- 
heated clouds. The dashed line represents the relative contribu- 
tion of directly-heated clouds to the torus emission at given dis- 
tance r from the AGN in the mid-plane, fun- The dash-dotted 
line shows the relative contribution (1 - /dh) of indirectly- 
heated clouds at r if fm = 1 . The solid line presents the actual 
contribution of indirectly-heated clouds due to the decreasing 
strength of the diffuse radiation field, /ih(1 - /dh)- 



case fill = P()(r,z) - 1 (and /dh = 0), the cloud's contribu- 
tion from indirect heating corresponds to the indirectly-heated 
database cloud emission at the same distance, tdb = fci, other- 
wise roB = r^i/ V/h- 

An interesting aspect is the actual contribution of indirectly- 
heated clouds to the overall torus emission. On one hand, with 
increasing r the fraction of clouds which are indirectly heated 
increases due to the increase in obscuration. On the other hand, 
the diffuse radiation field becomes weaker due to the absence of 
directly-heated clouds. In Fig. |2] we illustrate the relative con- 
tribution of directly- and indirectly-heated clouds to the torus 
emission at a given distance r in the torus mid-plane (i.e. z = 
or = 0). Here, we fix No = 5, a = -1.0, and R^ut = 5(Q. The 
relative contribution, fon, of directly-heated clouds at distance r 
from the AGN is shown as a dashed line in Fig.|2l As a dashed- 
dotted line, we show the relative contribution of indirectly- 
heated clouds if we assume that the vicinity of the indirectly- 
heated cloud is fully filled-up with directly-heated clouds to 
form the diffuse radiation field, i.e (1 - /dh)- This was assumed 
as the contribution of indirectly-heated clouds in iHonip et"al] 
(l2006l) . leading to the overestir nation of the SEDs at wavelengths 
> 20 fim (see Appendix A.2 in' Honig et alJl200^ . In reality, the 
number of directly-heated clouds around an indirectly-heated 
cloud decreases with r, as outlined in the previous paragraph. 
The solid line illustrates this effect. It shows the contribution of 
indirectly-heated clouds at distance r, considering the decrease 
of directly-heated clouds in the vicinity, /ih(1 - /dh)- As can 
be expected from energy conservation, this line is lower than 
the directly-heated cloud contribution and approaches it asymp- 



totically with increasing r. Given their general higher temper- 
atures, this illustrates that directly-heated clouds dominate the 
torus emission over indirectly-heated clouds - at least in more 
or less face-on geometries where line-of-sight obscuration ef- 
fects play a minor role. 

Finally, we want to justify why a single iteration of the 
diffuse radiation field is sufficient. Since the directly heated 
clouds are also indirectly heated, one would first have to sim- 
ulate directly-heated clouds. From these clouds, the diffuse radi- 
ation field has to be calculated and be included in the simulations 
of directly-heated clouds. Then, a new diffuse radiation field has 
to be determined and the process be repeated. After several it- 
eration steps, the final diffuse radiation field is obtained. The 
Monte Carlo method we used is quite flexible in implement- 
ing such a scheme without the need to reconsider direct illu- 
mination or previous steps (Kriigel 2008, p. 283): After the i-th 
simulation, a "corrected" diffuse radiation field has to be added 
which corresponds to the change with respect to the previous 



^v'drf- However, 



(;■ - l)-th iteration step, d^^.^^p - -^.drf 
it has be en noted that the diffuse radiation field co nverges quite 
fast (e.g. iHonig et al.ll2006; Nenkova et alj|2008ah . so that one 
can consider that the source function of the diffuse radiation 
field Sy-DRF 



.(0) 
' v;DRF 



£ Sy-DHi(l>)sin(pd(p, where Sy-ou is 



^ For a = -1.0, the normalized radial distribution function becomes 
I,. = (rlni?ou,r'. 



the source function of (AGN-only) directly-heated clouds, with- 
out the need of any further iteration. This makes simulations of 
cloud databases very time efficient. 



3. Emission from individual dust clouds 

In this section we will discuss dust cloud modeling results from 
step I of the torus simulation (see Fig.[T]). First, we discuss phys- 
ical properties of the dust clouds in the torus. Next, we detail 
the Monte Carlo simulation of individual dust clouds and ex- 
plain underlying assumptions. Finally we present cloud SEDs 
for both directly- and indirectly-heated dust clouds for different 
dust compositions. 



3.1. Properties of the dust clouds 

3.1.1. Observational and theoretical constraints on physical 
properties 

One fundamental property of the dust clouds in radiative trans- 
fer simulations is their optical thickness. If they are optically 
thin, i.e. ry <K 1, then they could be considered to have uni- 
form temperature and the thermally re-emitted flux can be ap- 
proximated by Fy - jiTyBy{T), where By(T) is the Planck func- 
tion at temperature T and Ty is the frequency-dependent optical 
depth. However, from observations of type 2 AGN, we know that 
both the dust and hydrogen columns are usually quite large, with 
the hydrogen column even Compton-thick. Associated torus dust 
temperatures are rather cool, pe aking atT ~ 200 - 300 K (e.g. 
iJaffe et al.ll2004t iTristram et al.ll2007) . Moreover, X-ray column 
density variability as seen in many AGN point to rather small 
clouds and column densities of the order of 10^^ cm"^ per cloud 
(e.g. Risaliti et al. 2007). To combine silicate absorption features 
as seen in type 2 AGN, Compton-thick obscuration, and small 
clouds, we infer that there are on average about 5-10 optically- 
thick clouds (tv » 1) along a line-of-sight to the center in a type 
2 AGN and that the volume filling factor of the torus (^y < 1, 
i.e. the torus is clumpy. These cloud properties are, at least qual- 
itatively, in agreement with models of self-gravitating dust and 
gas clouds in the shear of the gravitational potential of a super- 
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Table 2. Main physical properties of the dust compositions used 
in our study. 



Name identifier 


fractional composition 


^min 


^max 






silicates graphite 


dim) 


(pm) 


(pc) 


standard ISM 


0.53 0.47 


0.025 


0.25 


0.9 


ISM large grains 


0.53 0.47 


0.1 


1.0 


0.5 


Gr-dominated 


0.30 0.70 


0.05 


0.25 


0.7 



— Notes: " Sublimation radius rsub at a dust temperature Tsub = 1 500 K 
for Lboi = 10** erg/s. Sublimation radii for any other luminosity can 
be calculated by r = j-jub x (Lboi/10'** erg/s) Lboi is defined as 
^boi = f Lydv of the accretion disk. In our definition of the accre tion 
disk spectrum, the reverberation radius from , Kishimoto et ^ ( l2007h for 
Lboi = lO** erg/s would be r^j. = 0.38 pc, with approximately 25% scat- 
ter. 



massi ve black hole (iBeckert & Duschlll2004t iHonig & BeckertI 
l20Q7h . 

3.1.2. Dust composition 

Previous torus studies used dust properties in line with standard 
grain sizes and chemical composition (e.g. 'Scha rtmann et al.l 
[2005; Honig etal. 2006; Nenkovaetal. 2008b). In general, 
MRN distrib utions (referring to Mathis, Rumpl & Nordsieck; 
iMathis et al.l [l977) with a grain size power law distribution 
oc a~^'^ (a . . .grain size) was assumed with lower and upper 
limits flmin ~ 0.005 - 0.01 //m and amax ~ 0.25 /im, respec- 
tively. The chemical compositions invoked 47% graphite and 
53% s ilicates, either based on lDraind (l2003l) or lOssenkopf et alj 
(Il992l for silicate dust) optical properties. Such dust composi- 
tions are considered as "standard ISM". Recently it was sug- 
gested that the dust composition around AGN might be deviat- 
ing from standard ISM. While Suganuma et al. (2006) confirmed 
the L'^^-dependence of the sublimation radius, the measured K- 
band reverberation radii of type 1 AGN are a fa ctor of 3 smaller 
than w hat is expected from standard ISM dust. iKishimoto et alJ 
ilOO 'f) suggest that this discrepancy might be caused by domi- 
nation of larger grains, at least in the inner part of the tori. 

We aim for testing the impact of different dust compositions 
on the cloud SEDs. We synthesize standard ISM mixes invok- 
ing 47% graphite (using opt ical constants from Draine (2003)) 
and 53% silicates using both lDraili3 (12003) and Ossenkopf et al. 
( 11992) optical constants. In order to explore recent suggestions, 
we also create a dust mix with standard ISM grain size distri- 
bution but only containing grains between 0. 1 yum and 1 jim in 
size ("ISM large grains"). A last composition explores SEDs of 
dust which is dominated by intermediate to larger graphite grains 
(70% graphite, 30% silicates; "Gr-dominated"). The composi- 
tions and their respective properties are summarized in Table|2l 

3.2. Monte Carlo simulations of dust clouds 

Since the exact dust distribution within a cloud is not known, 
we use a uniform density throughout the cloud. While this might 
be oversimplified, we note that independent of the actual density 
distribution of an optically thick cloud, the bulk of the cloud re- 
emission originates from the optically-thin surface layer of the 
hot, directly-illuminated face of the cloud. Moreover, the gradi- 
ent of the temperature at distance r from the cloud surface de- 
pends mostly on the total optical thickness at r and is supposedly 
very similar in all density distribution laws once the same optical 



depth is reached. So, the temperature of the non-illuminated side 
of the cloud is more or less independent of the dust density law 
within the cloud for the same optical depth through the diameter 
of the cloud (="total optical depth" Ty). 

The cloud is subdivided into grid cells, where the dust den- 
sity within each cell is optically thin. This guarantees that the 
steep temperature gradient in the surface layer on the hot face of 
the cloud is properly sampled. 

For our cloud simulation, we use a spherical cloud geometry. 
In principle any kind of shape could have been chosen. However, 
the main difference between a sphere and, say, an ellipse is the 
different ratio of radial-to-perpendicular optical depth. As men- 
tioned above, if the cloud has a total Ty » 1, then the dominant 
emission is coming from the optically-thin layer on the hot side 
of the cloud. Since the temperature in the layer is mostly inde- 
pendent of the density distribution in the rest of the cloud, the 
overall temperature of the hot cloud face is similar. For the cold 
face of the cloud, the main criterion is again the total optical 
depth within the cloud. For comparable total optical depths of 
two clouds with arbitrary shapeQ, the cold-side temperatures are 
the same and, thus, the emerging emission. In summary, unless 
one insists on very peculiar cloud surfaces, the exact shape does 
not matter too much and a spherical cloud presumably catches 
the essence of the dust clouds in AGN tori. 

For time-efficient simulations, we optimize the Monte Carlo 
code for the cloud-AGN-configuration. For that, we first de- 
termine the sublimation radius r^ub - fiTsub) for a given dust 
chemistry and grain size composition, assuming a dust subli- 
mation temperature T^ub - 1 500 K. Since rsub scales with the 
AGN luminosity L as rsub <^ i'^^, the incident AGN flux on 
a dust cloud at rsub is independent of the actual luminosity, 
■P'AGN('"sub) ^/'"iib ~ constant. Thus, a single cloud database 
can be used for all AGN luminosities. Since the dust clouds are 
believed to be small, with a radius-to-distance-ratio Rc\/r ^ 1, 
we assume that the incident AGN radiation (see Sect. |A|l enters 
the model sphere on parallel rays. Moreover, the rotational sym- 
metry of our clouds allows to reduce the model space for Monte 
Carlo simulations to a planar 3-dimensional grid with only 1 grid 
cell in z-direction. These two modification allow for a maximum 
number of photon packages to be simulated in a short time with- 
out suffering from strong Monte Carlo noise in the temperature 
of the individual grid cells. To be even less dependent on Monte 
Carlo noise, the final cloud SEDs and images are calculated via 
raytracing of the emission from each cell into the observer's di- 
rection (i.e. for different phase angles of the cloud): First, the 
temperature distribution is reconstructed from the temperature 
grid; then, the source function (thermal re-emission and scatter- 
ing) for each cell is calculated and the path through the cloud 
is followed out of the cloucfl In this way the angle-dependent 
source function of the cloud S y{(p) is time-efficiently recovered. 
The described simulation method is the same for both directly- 
and indirectly-heated clouds. The only difference occurs in han- 
dling the incident radiation: in the case of directly-heated clouds, 
the AGN spectrum is used while for indirectly-heated cloudes 
the diffuse radiation field is the heating source. 



^ The actual total r^'s that need to be compared depend on the exact 
shapes/elongations. Thus, one might see the total Ty of the cloud in our 
model as some mean value representing all kind of shapes. 

* Note that this treatment limits the handling of scattering to single 
events. For our optically-thick clouds we tested this method against rig- 
orous treatment of scattering with the mcsim code and found that any 
noticeable impact of the single-event-assumption is limited to wave- 
lengths i 1 ^m and is of the order of 2%. 
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ISM standard (D03) ISM standard (Oss.) ISM large grains Gr~dominated 




0.1 1 10 100 0.1 1 10 100 0.1 1 10 100 0.1 1 10 100 

wavelength (/xm) wavelength (/i.m) wavelength (/^m) wavelength {fim) 



Fig. 3. Model source functions 5 v of directly AGN-heated clouds for several different dust configurations (from left to right column: 
ISM standard configuration with Draine (2003) silicates, Ossenkopf et al. (1994) silicates, ISM large grains, and Graphite-dominated 
dust). Each row shows the emission of a cloud with Ty - 50 at the given AGN distance r = rsub (top), 10 r^ub (middle) and 100 rsub 
(bottom), scaled to the total incoming flux of the AGN, Fagn - iboi/(47rr^). The corresponding hot-side temperature in units of 
Tsub - 1 500 K are also given. The various lines represent a specific phase angle, (p, of the cloud in steps of 30° (red: = 0°, i.e. full 
hot side; purple: = 180°, i.e. full cold side). For comparison, we overplotted the corresponding black-body emission By(T) with 
the same temperature as the hot side of the cloud at the same distance from the AGN (dashed line). 



In summary, our simulations of dust clouds pr o ceeds as 
follows (see also iHonig et alJ (l2006t) : iNenkova et aP (l2008ah ): 

(1) Monte Carlo simulations of directly AGN-heated clouds; 

(2) Approximation of the diffuse radiation field (i.e. heating of 
clouds not directly exposed to the AGN) by averaging the lo- 
cal emission of directly-heated clouds (see also Sect. 12.31 for a 
description of the assessment of the local circumstances within 
the torus); (3) Monte Carlo simulation of clouds heated by 
the averaged diffuse radiation field. This results in two cloud 
databases which are used for the torus simulations: one for 
directly-illuminated clouds and one for indirectly-heated clouds. 



3.3. Resulting cloud SEDs 

In the following we present model source functions of directly- 
and indirectly-heated clouds. Simulations have been carried out 
for the standard-ISM, ISM large-grains, and Gr-dominated dust 
configurations as introduced in Sect. 13.1.21 



3.3.1. Directly-heated clouds 



Standard ISM The starting point for our dust composition study 
will be the "standard ISM". Physical dust properties are listed in 
Tab. |2] We show two versions of the standard ISM configura- 
tion: one u ses Draine (2003) silicates, the other one consists of 
lOssenkopf et al.l (11992) silicates (see Sect. 13.1. 2] |. The most im- 
portant difference between both silicate types concerns the cen- 
tral wavelength and width of the silicate feature at around 10 /zm. 
For a given size distribution, the Draine silicate feature is much 
wider than the Ossenkopf feature. In Draine silicates, the central 
silicate feature wavelength is at ~ 9.5 ;um while the Ossenkopf 
silicate feature is centered at 10.0 fim. The 10.0 fim central wave- 
length seems to be more consistent with peak silicate emission 
features of type 1 AGN as observed with Spitzer, while absorp- 
tion features in type 2s usually have their center at 9.7 fim (e.g. 
IShi e t al. 2006). We note, however, that radiative transfer effects 
can affect the exact shape of the feature (but see discussions in 
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Fig. 4. Wavelength-dependent emission ratio between the source function of indirectly-heated clouds S^xw and the hot face of 
directly-heated clouds 5v.;dh(<^ = 0°). The various lines represent different distances r from the AGN with r - \, 1.5, 2.5, 4, 7, 10, 
15, 25, 40, 70, 100 and 150rsub from top-red to bottom-blue. The left panel shows simulation results for "ISM standard" dust with 
Ossenkopf et al. silicates, the middle panel presents "ISM large grains", and the right panels shows "Gr-dominated" dust. 



iNikutta et al.ll2009t iHonig etani20Tol: iLandt et al.ll20ia which 
disagree on the details of this effect). 

Simulation results for directly-illuminated clouds with stan- 
dard ISM dust configurations are shown in the two left columns 
of Fig. [3] We present clouds with a total optical depth ^^^ - 50 
at three different distances r - rjub, 10 x rjub, and 100 x rsub 
from the AGN. Along with the distances, we also provide the 
maximum temperature of the hot side of the cloud in units of 
the sublimation temperature Tsub - 1 500 K. The overplotted 
dashed lines show the black-body source function By{T) with 
the same temperature as the hot side of the clouds at the same 
distance from the AGN. For both types of silicates, the temper- 
ature gradient on the directly-illuminated face is similar, with 
some trend of slightly higher temperatures for the Ossenkopf 
dust. The r = 10 x rsub clouds serve best for illustration since the 
near-IR color is slightly redder for Draine silicates. The different 
temperature gradients are a result of slightly different absorption 
efficiencies - or more precisely: Planck mean opacities which, 
at the end, determine the temperature at a given radius in the ra- 
diative transfer equations. The cloud emission is lower than the 
corresponding black-body emission by more than a magnitude 
indicating that ISM standard mixture clouds are not well rep- 
resented by any black-body approximation. This is true for all 
dust compositions discussed here. It is worth noting that for the 
standard ISM configuration, the presented clouds with ry = 50 
are mostly optically thin in the infrared. This leads to only small 
differences in SED shapes for different phase angles (/0 so that 
the hot-side emission shines all through the cloud and signifi- 
cantly contributes to the cold-side emission, e.g. as seen by the 
9.7/10/vm silicate emission feature in the (f> = 180° SEDs. 

ISM large grains While near-IR reverberation mapping results 
of type 1 AGN confirmed the L'^^ -dependence of rsub, the ob- 
served /T-band reverberation mapping radii r^^, are smaller by 
approximately a factor of 2-3 than what has been inferre d 
from ISM dust dSuganuma et all 120061: iKishimoto et al.ll2007h . 



^ The phase angle <p is defined as the AGN-cloud-observer angle. <p = 
0° denotes the full visibility of the directly-illuminated face of the cloud 
("full moon"), = 180° describes the situation when only the cold side 
of the cloud is seen ("new moon"). 



Suggestions to solve this problem involve larger sizes of the 
dust grains than used in standard ISM distributions - at least 
in the innermost part of the torus which dominates the IR emis- 
sion of type 1 AGN. To test this suggestion, we simulated dust 
clouds for a standard ISM composition with Ossenkopf sili- 
cates and limited the grain sizes to 0. 1 - 1 jum. This "ISM large 
grain" composition has a sublimation radius rsub = 0.5 pc x 
(Lboi/lO^*' erg/s)'^^ which is closer to the observed r^^, than the 
standard ISM dust (see Tab.|2]i. 

The simulated dust cloud SEDs for the ISM large grain com- 
position are shown in the third column of Fig. |3] The hot-side 
temperature gradient is comparable to clouds with ISM standard 
dust. The most notable difference concerns the change of the 
SED shape from (p - 0° to 180°. While the previously-explored 
clouds show quite similar SEDs for all phase angles, the ISM 
large grain clouds become very red at near- and mid-IR wave- 
bands for increasing cp. Moreover, the silicate feature at IQfim 
turns from emission into absorption. The reason for this differ- 
ence is the optical depth in the infrared: For both silicate and 
graphite, large grains have a much flatter extinction efficiency 
in the infrared. Thus, the cloud-internal extinction in the IR in- 
creases. While Ti.i/jm/Ty ~ 0.05 - 0.06 for the ISM standard 
configuration (depending on the silicate species), the ISM large 
grains have T2.2fim/Tv ~ 0.9. In the silicate feature at 9.1 /lOj-im, 
the ratios are rsi/ry ~ 0.04 and 0.1, respectively. As a conse- 
quence, a cloud with ry - 50 is also optically thick in the near- 
and mid-infrared. Thus, when seeing the cloud under large (p, the 
emission from the hot side of the cloud is blocked and the emerg- 
ing SED is dominated by cool dust emission. This will also have 
an effect for the simulated torus SED since IR radiation can be 
efficiently absorbed which presumably leads to redder IR colors 
(see Sect. 14. 1.1b . 

We note that the difference between black-body and hot-side 
cloud emission is smaller than for the other dust composition, 
as seen in Fig. [3] This implies that large dust grains are better 
represented by a black-body approximation than smaller grains, 
although the difference is still almost an order of magnitude. 

Graphite-dominated dust In general, silicate dust grains have 
lower sublimation tem peratures than graphite grains (e.g. 
ISchartmann et alJ I2005L for simulations with AGN radiation). 
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This might resuh in a dearth of silicates in the torus, at least 
in the inner part. Moreover, since AGN are rather strong X-ray 
emission sources, it is doubtful that smallest dust grains can eas- 
ily survive without being photo-destructed. For that, we modi- 
fied the ISM composition by using only grains > 0.05 /urn and 
changing the chemical mixture to 70% graphite grains and 30% 
Ossenkopf silicates. 

In the right-most column of Fig. [3] we show simulated cloud 
SEDs for this "Gr-dominated" dust. The overall SED shape re- 
sembles those of the standard ISM configurations. The dust is 
optically thin in the IR with T2.2nm/Tv - 0.06 and rsi/ry - 0.02. 
However, two important differences can be seen. First, the sili- 
cate emission feature at lOyum is significantly reduced due to 
the reduction of silicate dust. This might be interesting given 
the relative weakness of silicate absorption and emission fea- 
tures observed in Seyfert galaxies (but see Sect. 14. l.T] for actual 
simulations). Second, the temperature gradient outward from 
rsub is slightly flatter. As can be seen in the middle panel, the 
near-IR emission in the Gr-dominated dust clouds is bluer than 
for the standard ISM configurations and the ISM large grains. 
Therefore, it can be expected that for a given set of torus pa- 
rameters (see Sect. 12.21 ). the radial emission size in the near- and 
mid-IR is slightly different than for the other dust compositions. 



3.3.2. Indirectly-heated clouds 

In Fig. m we illustrate the contribution of indirectly-heated 
clouds. The left panel shows the ratio of indirectly-heated cloud 
emission to directly heated emission as the ratio of the source 
functions S v-jh/S v^DH{<f> = 0°) for the ISM standard dust with 
Ossenkopf silicates at different distances r from the AGN. In the 
middle panel, we present results for the ISM large grain con- 
figuration, and in the right panel Gr-dominated dust was used. 
Except for the clouds closest to the AGN (r ^ 5 x rsub), the emis- 
sion contribution of indirectly-heated clouds in the near- and 
mid-infrared below lOjum is <K 0.1 for the ISM standard and 
Gr-dominated dust. This means that indirectly-heated clouds are 
much cooler than the corresponding directly-heated cloud at the 
same r. This is slightly different for the ISM large grain configu- 
ration: Since these clouds are also optically thick in the infrared, 
they absorb the incident diffuse radiation field - which is domi- 
nated by IR photons - much more efficiently than the other two 
dust configurations. This results in better indirect heating so that 
the temperatures become higher 



4. Torus model results and discussion 

In this section we will present results from the full torus sim- 
ulations as described in Sect. 12.1! and illustrated in Fig. [T] 
There have already been a number of SED parameter studies 
on clumpy torus rnodels i n hterature jP ullemond & van Bemmej 
2005- 'H onig et al.ll2006l:ISchartmann et al.n200 8: Nenko va et all 
2008b), so that we want to focus on the dust distribution and how 
it connects to observations. In particular, we will show that cer- 
tain SED and interferometric properties are almost exclusively 
depending on the radial distribution of the dust clouds in the 
torus. The 0-distribution will be used to discuss parameter de- 
pendencies, but the z-distribution will be included in the conclu- 
sions drawn from the SED studies in Sect. 14.1.51 Both distribu- 
tions yield very similar results. 



4.1. Torus Spectral Energy Distributions 

All torus SEDs will be scaled independent of AGN luminosity 
and spatial scaling. The torus model SEDs vFy shown in this 
section represent the emission at a distance corresponding to 
the sublimation radius rsub of the respective dust composition. 
Observed fluxes vfy can be easily compared to model fluxes vFy 
via the relation 

vFy = vfy ■ i^] (7) 

\'"sub/ 

where Dl is the luminosity distance to the AGN and r^ub is the 
sublimation (or reverberation) radius for the modeled dust com- 
position which can be taken from Table|2] By that, a comparison 
of observed and model SEDs allows for constraining the subli- 
mation radius from which, if needed, the bolometric luminosity 
can be inferred (see Table|2|. 

4.1 .1 . The radial distribution of dust clouds 

In the parametrization of the torus the radial distribution of the 
clouds takes a key role in determining the emission, as illus- 
trated in eqs. (O, Q, and This is reflected by the differences 
in simulated SEDs when varying the power-law index a in the 
radial dust-cloud-density distribution //, r". In Figs.|5]&|6]we 
present comparison of model SEDs for different radial power- 
law indices a = 0.0, -0.5, -1.0, -1.5, and -2.0, indicated by 
colored lines in the plot. Fig.|5]represents a face-on line-of-sight 
onto the torus as in type 1 AGN (/ = 0° while Fig.|5]shows SEDs 
for a type 2 case (i.e. face-on; / = 90°). Each row represents a 
different dust composition as discussed in Sect. 13.1! (top row: 
Ossenkopf et al. standard ISM; middle row: ISM large grains; 
bottom row: Gr-dominated dust). Each column shows variations 
of the mean number of clouds in the equatorial line-of-sight 
A^o (see Sect. 12.21 ). increasing from left to right (A^o = 2.5, 5, 
7.5, and 10). Other fixed parameters are 7?out = 150, 6q = 45°, 
R,i,o = 0.035, = 1, and = 50. 

All models have in common that the total SED becomes red- 
der when the power-law becomes flatter, simply because cooler 
dust at larger distances is involved with the flat power-law distri- 
butions. The standard ISM and Gr-dominated dust distributions 
show similar overall SEDs and similar trends when changing A^o 
or a. For small A^o, the silicate emission features are strongly 
pronounced, and there is quite some difference in the continuum 
shape for different a. While the flat distributions (a = 0.0 to 
-1.0) peak in the continuum at long mid-IR wavelengths, the 
steeper distributions (a - -1.5 and-2.0) have their continuum 
emission peak in the near-IR, which is caused by the fact that 
in the latter case most of the dust is confined to small distances 
from the AGN, i.e. the average dust temperature is rather high. 
This is the same trend for type 1 and type 2 line-of-sights, al- 
though type 2 SEDs appear slightly redder than their type 1 
counterparts. 

A special case in this study are SEDs for the ISM large grain 
dust. While they follow the general trend of redder SEDs for 
shallower dust distributions, the peak continuum emission is lo- 
cated at wavelengths > 10 //m in vFy for all selections of a and 
A^o in the type 2 case, and their continuum is generally redder 
in type 1 orientations. The reason for this behavior is the dust- 
specific T in the infrared (see Sect. 13.3.1b . For the clouds with 
total Ty - 50, the large grains are optically thick in the near- 
and mid-IR (tir > 1), while the other dust compositions are 
optically thin in the mid-IR. This can also lead to stronger statis- 
tical effects, depending on the actual distribution of clouds (see 
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Fig. 5. Model SEDs of a face-on type-1 AGN (inclination angle ; = 0°). The top row shows models with the standard ISM compo- 
sition, the middle row is for large grains, and the bottom row represents GR-dominated dust. From left to right, the columns show 
an increasing mean number of clouds in the equatorial line-of-sight, A^o = 2.5, 5, 7.5, and 10, respectively. In each panel, we plot 
SEDs of one random cloud distribution for radial power law indices a = 0.5 (red), -0.5 (orange), -1.0 (green), -1.5 (light blue), 
and -2.0 (dark blue). Fixed parameters are Roui = 150, 6q - 45°, Rd^u - 0.035, b = I, and Td = 50. 



iHonig et"aDl2006h . Considering overall similarities in the SEDs 
and in the silicate features of most type 1 cases, it seems difficult, 
however, to constrain the chemistry by SED observations. 

While the overall SED shows considerable dependence on a, 
it does not change significantly when varying A^o for both type 1 
and type 2 line-of-sights. There is only a small change in color 
from A'o = 2.5 to 10 which can be noted as a smaller dispersion 
between the extreme a-cases by comparing the left-most with the 
right-most panels in Figs.|5]&|6] In fact, the red continuum in flat 
distributions becomes bluer and the blue, steep dust distribution 
become redder This is caused by obscuration effects within the 
torus: if A^o increases, clouds at larger distances contribute less 
to the overall emission since there is less chance that they are 
directly exposed to AGN emission. That effect is strongest for 
the flat distributions. In total, however, torus-internal obscura- 
tion effects on the SED shape, originating from the selection of 
No, are minor. This can be taken as a justification of the sim- 
plified clumpy torus model for type 1 AGN as introduced by 
[Kishimoto et al. (2008) and quantified in eq. (|6]l which captures 
the essence of the (overall) brightness distribution in pole-on ge- 
ometries. However, the selection of A^o seems to have an impact 



on the strength of silicate feature, which will be discussed in the 
next section. 

4.1.2. Average number of clouds along an equatorial 
line-of-sight 

In the previous Sect. 14.1.11 it has been shown that the effect of 
varying A^o on the overall SED is minor. On the other hand, the 
silicate feature in both type 1 and type 2 AGN changes from 
strong emission to weak emission or absorption with A^o increas- 
ing from 2.5 to 10. 

In Figs.|2]&[8] we show the lOjum silicate features for differ- 
ent model parameters for type 1 and type 2 AGN respectively. In 
the case of type 1 AGN we show feature strengths for Ossenkopf 
ISM standard, large grain, and Gr-dominated dust mixes. To iso- 
late the features, we followed t he continuum spline fit proce- 
dure proposed bv ISirockv et alj (2008). For that a cubic spline 
has been fitted to the regions at 5-7 jum, 14-14. 5jum, and 25- 
31.5 yum. The fitting result is considered representative for the 
"continuum" underlying the silicate features at 10 and 18//m. 
While this method is applicable to Spitzer data, ground-based 
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Fig. 6. Model SEDs of an edge-on type-2 AGN (inclination angle / - 90°). The top row shows model models with the standard ISM 
composition, the middle row is for large grains, and the bottom row represents GR-dominated dust. From left to right, the columns 
show an increasing mean number of clouds in the equatorial line-of-sight, A^o = 2.5, 5, 7.5, and 10, respectively. In each panel, we 
plot SEDs of one random cloud distribution for radial power law indices a = 0.5 (red), -0.5 (orange), -1.0 (green), -1.5 (light 
blue), and -2.0 (dark blue). Fixed parameters are Rout = 150, 6q - 45°, Rci;[) - 0.035, b = 1, and t^i = 50. 



data usually have a much s maller wavelength c overage due to 
the atmospheric cutoffs (e.g. iHonig et al.ir20ld ). For such data 
we recommend the use of a linear fit to the 8.5 and 12.5 yum 
fluxes as continnum representation. These wavelengths should 
be reasonably unaffected by the lOyum silicate feature - at least 
at typic ally obser ved strengths (e.g . Mason et al. 2006, 2009| 
iHorst et al. 2008; Honig et al.ll2010l) . As in the previous figures, 
we present increasing A^o from 2.5 to 10 in the columns from 
left to right. In each column, we color-coded the individual a- 
indices from -0.0 to -2.0 as in Fig.|5] From top to bottom, the 
rows show standard ISM, large grain, and Gr-dominated dust 
compositions. In general, the ISM standard configuration shows 
stronger emission features than the other dust configurations. 

Most of the change in the strength of the silicate feature hap- 
pens when varying A^o- Higher A^o corresponds to weaker emis- 
sion features in type Is. The same applies to type 2 AGN, how- 
ever with a tendency of weaker emission features and slightly 
more pronounced absorption features than in type Is. In addi- 
tion to the influence of A^o on the feature strength, shallower 
dust distributions show less pronounced silicate emission fea- 
tures and stronger absorption features than steeper dust distri- 



butions. This is due to the interaction of obscuration in vertical 
direction (for larger A^o, more clouds are present in vertical direc- 
tion) and cloud distribution (flatter distributions are dominated 
by cooler dust emission). 

In summary, while the overall SED depends only on one pa- 
rameter (for a given dust distribution; as shown in the previous 
section), the exact strength of lOjjm silicate feature is depending 
on at least A^o and a, with some emphasis on the former. Other 
parameters which may have an influence on the silicate feature 
are discussed in the next sections. 

The standard ISM dust configuration with Td around 50 
seems to capture the essence of a typical torus SED including 
moderate silicate emission in type 1 and absorption in type 2 
line-of-sights for a broad range of model parameters. We will, 
thus, use the ISM standard dust with Ossenkopf silicates in the 
following sections. 

4.1 .3. The opening angle of the torus 

A potential source of uncertainty in determining A^o from the sil- 
icate feature (and maybe a from the continuum SED) could be 
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Fig. 7. Continuum-normalized lOyum model silicate features in type-1 AGN (inclination angle i - 15°). The top row show model 
models with the standard ISM composition, the middle row is for large grains, and the bottom row represents Gr-dominated dust. 
From left to right, the columns show an increasing mean number of clouds in the equatorial line-of-sight, A^o = 2.5, 5, 7.5, and 10, 
respectively. In each panel, we plot SEDs of a representative random cloud distribution for radial power law indices a = -0.5 (red), 
-1.0 (yellow), -1.5 (green), and -2.0 (blue). Fixed parameters are R^^a = 150, 6*0 = 45°, Rd-o = 0.035, b ^ I, and t^i = 50. 



the geometrical thickness of the torus that we parametrized as the 
half-opening angle Oq. In Fig.|9]we show how changing Oq affects 
the mid-lR spectral index a and the silicate feature for standard 
ISM dust. For that we plotted the "apparent" optical depth rsi in 
the silicate feature as determined in the previous section versus 
the spectral index a. The mid-IR spectral index a has been deter- 
mined from the continuum spline fit as described in Sect. 14.1.21 
We then made a linear fit to the 7.0-8.5//m and 13.9-14.6yum 
region to determine the spectral slope in the A^-band from which 
a has been calculated. Within Fig. |9] we show the location of 
models ranging from a = 0.0 to -2.0 and with A^o - 2.5 - 10 for 
00 - 60°, 45°, and 30°. The left panel shows results for a pole-on 
view on the torus, the right panel represents a type-2-like edge- 
on view. Varying Oq mainly affects the strength of the silicate 
feature in a way that emission features are weakened when the 
half-opening angle is smaller. The reason is that smaller Oq re- 
sults in more clouds along vertical line-of-sights which on one 
hand obscure the emission features from the torus mid-plane in 
type 1 AGN and on the other hand preferentially expose their 
cooler sides (= less strong emission features) to the observer. 
This effect is less pronounced for type 2 AGN where clouds af- 



fected by obscuration are those on the far-side of the torus which 
would be seen almost with their full hot emission side. The over- 
all a and Tsi varies less in type 2 line-of-sights when varying Oq. 

In summary, varying 0q alters primarily the silicate fea- 
ture depth in type 1 cases while the spectral index is al- 
most unaffected. Thus, using observables a and rsi leaves 
us with a degeneracy in model parameters A^o and 9q. 
However, sample studies showing type-l /type-2 rat i os of 
1:1 to 1:3 (e.g. iM aiolino & Rieke 199^ iLacv et all 12004 
lMartmez-Sansip"e et al. 2006) and observed narrow-line region 
opening cones around 90° in many nearby objects favor Oq ~ 
40-50°. For indiv idual objects, Oq may be smaller or larger (e.g. 
iMason et"al ] 12009) than this average. When considering, how- 
ever, a larger sample of AGN, the use of 6*0 = 45° as the sample 
average is well justified to break the A^o -6*0 -degeneracy and get 
an idea what range of A^o is covered by the AGN in the sample. 
On the other hand, the uncertainty in the individual clouds' op- 
tical depths remains in such an analysis unless independent ev- 
idence narrows down t^i (e.g. general agreement of SEDs with 
unification scheme, IR visibilities. Hydrogen column densities, 
theoretical considerations). Thus it seems more viable to give 
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Fig. 8. Continuum-normalized 10/im model silicate features in type-2 AGN (inclination angle / = 90°) using the ISM standard dust 
composition. From left to right, the panels show an increasing mean number of clouds in the equatorial line-of-sight, A^o = 2.5, 5, 
7.5, and 10, respectively. In each panel, we plot SEDs of a representative random cloud distribution for radial power law indices 
a = -0.5 (red), -1.0 (yellow), -1.5 (green), and -2.0 (blue). Fixed parameters are Rom = 150, do = 45°, Rao = 0.035, b = I, and 
Tcl = 50. 
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Fig. 9. Dependence of observed silicate-feature optical-depth tsi and mid-IR spectral index a on model parameters. The grids show 
how Tsi and a change with a (dotted lines; from a - 0.0 to -2.0) and A^o (dashed lines; from A^o = 2.5 to 10) for different half- 
opening angles of the torus (red: 60°, green: 45°, blue: 30°). Left: Type 1 AGN with / = 0°; Right: Type 2 AGN with / = 90°. The 
grid lines are results for one random arrangement of clouds for the given set of parameters. 



a sample A^o-T"ci-range or providing the A^o constraint under the 
assumption of a certain Td. 

4.1 .4. The outer radius of the torus 

When modeling the IR SEDs, there is often a question arising: 
How far outward does the torus extend in the model? Since our 
model defines Rg^ as a parameter, one might think that the /?out 
value would be the correct answer However, it is not all that 
simple. Due to the decrease of cloud temperature with distance 
and the corresponding decrease of total intensity, clouds at dif- 
ferent distances from the AGN contribute different fractions to 
the total torus flux at a given wavelength. In general, clouds at 



the innermost torus region contribute most of the near-IR light 
while clouds at larger radii are the dominant source of the mid- 
IR emission. Being more quantitative we can approximate 

(rri \ —2...— 2.8 

for dust grains in radiative equilibrium (e.g. iBarvaini ^ Il987h . 
where T^uh is the dust sublimation temperature (~1500K) and 
T is the dust temperature at distance r from the AGN. The range 
given for the exponent covers the range from black body grains 
(large, graphite) to typical ISM dust (see also Barvainis 1987). 
Based on this equation and Wien's law, it is possible to estimate 
the maximum radius which contributes to the torus emission at 
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a given wavelength. Thus, the maximum size of the emission is 
depending on the observed wavelength. On the other hand, the 
exact shape of the SED (and also the observed size of the emis- 
sion region; see Sect. 14.2b depends on the dust distribution (see 
previous sections), so that eq. ([Sjl provides only a rough esti- 
mate. However, choosing /?out smaller than the size correspond- 
ing to the temperature {- wavelength) of interest may result in 
an artificial cut-off in the brightness distribution of the models. 
For example, when selecting Rout - 25, the cut-off occurs at 
temperatures of T ~ 475 K for ISM dust which corresponds to 
wavelengths around 8 fim. 

The torus brightness distribution does not only depend on the 
dust temperature (or source function of clouds) but also on the 
dust distribution as shown in Sect. 12.21 Steep dust distributions 
have most of their dust at small radii so that they are mostly not 
sensitive to 7?out- However, the actual choice of Rom will become 
important for the mid-lR SED once the dust distribution is very 
flat (a ^ -1) or even inverted (a ^ 0)0. In this case the inte- 
gration of the brightness profile will not converge for Roui — » oo 
and the result depends on /?out (see eq. It is, thus, important 
to select Rom properly to avoid artificial cut-offs in intermedi- 
ate cases. Moreover for quite extended distributions, such as the 
presented a - 0.0 case. Rout has to be considered as a full model 
parameter. 

Recent IR interferometric observations and /iT-band rever- 
beration mapping measurements suggest that the torus bright- 
ness distributions are qu it e extended at least up to ~ 100 x r^ub 
dKishimoto et al."2009a''^, 'Honig et al."2010) in contrast to the 
discussion in Nenkova et al. (2008b). For a better constraint on 
the "real" Rout the use of ALMA at even longer wavelengths will 
be crucial. However, it may be difficult to set a strict limit of the 
torus (= AGN-heated dusty region) with respect to its galactic 
environement (= star or starformation-heated). 



4.1 .5. Using SEDs to determine tine radial power-law index a 

Summarizing the subsection on torus SEDs, we showed that the 
slope in the mid-IR, represented by a, is strongly related to the 
radial dust distribution in the torus. This seems to be the case 
despite several other parameters which potentially influence the 
torus emission. It is, however, consistent with expectations based 
on theoretical considerations (see Sect. 12.2b . 

In Fig. [To] we present a more quantitative approach to the 
relation between a and a in type 1 AGN. For each model a pa- 
rameter, we varied 0o (different symbols) and A^o (color-coded), 
and plotted the resulting a versus a. In addition we assess the 
influence of torus inclination in type 1 AGN by varying / from 
0° to 30° (offsetted in the plot). To be even more general we also 
included simulations for both the spherical (bluish colors) and 
cylindrical (reddish colors) vertical distributions (see Sect. 12.21 ). 
This results in 72 modeled a values for each a. The black filled 
circles with errors note the mean of these models for each a and 
the standard deviation around the mean valu^ These mean val- 
ues have been fitted by a straight line with slope -0.93, and the 
overall correlation of a vs. a has a Spearman rank of -0.97. The 
main outliers from the fit are geometrically very thick (0o - 30°), 
higher inclination (/ - 30°) cases with high A^o using the spher- 



^ a = - 1 seems to be a borderline case. Although the integration of 
the brightness profile does not formally converge, obscuration effects 
will reduce the dependency on R^ut if properly selected. 

^ For convenience we calculated the standard deviation from all 
model a values for each a value, although the distribution of a's around 
the mean is non-Gaussian, non- symmetric. 
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Fig. 10. Relation between the mid-IR spectral index a and the ra- 
dial dust distribution power-law index a in type 1 AGN. Model 
points are shown for all combinations of A^o = 2.5, 5, 7.5, and 10 
(color-coded symbols) and half-opening angle 30° (squares), 45° 
(circles), and 60° (diamonds) for both spherical (bluish-colored 
symbols) and cylindrical (reddish-colored symbols) vertical dis- 
tributions. In order to cover a range of type-1 inclinations, results 
for / = 0° (face-on; upward offset), i - 15° (nominal a value) and 
/ = 30° (downward offset) are presented for all a-values. The 
black-filled circles with errors denote the nominal a-value for 
each a and one standard deviation. The nominal fit using all val- 
ues (solid black line) has a slope of -0.93 and a Spearman rank 
correlation of -0.97. Fits using only results with half-opening 
angle 30° (dash-dotted line), 45° (dashed line; almost coincident 
with the solid line), and 60° (dotted line) are also shown for il- 
lustration. 



ical distribution, while the general relation is quite tight. In fact 
the slope of the fit does not change too much when fixing 6^ (see 
lines for 6*0 - 30° (dashe-dotted), 45° (dashed; almost coinci- 
dental with solid line), and 60° (dotted)), ranging approximately 
from -0.87 to -0.95. The nominal fit relation for the mean val- 
ues including errors is 



(-2.16 + 0.16)- (0.93 ± 0.1 n-ff 



(9) 



Within the range of errors illustrated in Fig.[TOl this relation may 
be used to extract information about the dust distribution from 
mid-IR spectra. 

4.2. Interferometry of dust tori 

Due to the mas-spatial scales involved in torus observations, it is 
difficult to spatially resolve the torus. The most promising tool 
for this task is IR interferometry, which w as already success- 
fully used for a r iumber of nearby objects jSwain et al] 2003 1; 
Wittkowski et al.' '2004; 'Weigeh et al.' 2004; Jaff'e et al. 2004[_ 
Poncelet et al. 2006; Meisenheimer et al. 2007; Tris tram et al 
2007; Beckertetal. 2008; 
I2009t iBurtscher et al.i,2009 



Raban et al. 2009; Tr istram et al 
,Kishimoto et al.ii2009ta) . However, 
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except for reconstructe d images of NGC 1 068 obtained from 
Speckle interferometry (IWeigelt et al.l 1200 4*). long-baseline in- 
terferometry usually provides spatial information coded as 
Fourier space observables, i.e. phase and visibility information 
for a given (projected) baseline and position angle. While the 
object's phase is important for image reconstruction, the visibil- 
ity contains the fundamental spatial information about the object 
(size, elongation). This information can be used as complemen- 
tary to total flux SEDs, in particular if for a given wavelength 
both total and correlated fluxes are available. 

The figures in the following subsections are simulated for an 
AGN at 15Mpc distance with an inner radius rsub = 0.025 pc. 
They can be used as predictions for any AGN by re-scaling the 
jc-axis baseline scale according to 

(Da-is] 

BLaGN - BL^odel ■ \ — —] (10) 
\ '"0.025 / 

where D^;ig is the AGN's angular-diameter distance in units of 
15 Mpc and ro 025 is the reverberation/sublimation radius in units 
of 0.025 pc. 

4.2.1 . The baseline-dependence of the visibility 

As shown in literature, adding interferometric information is 
very constraining for modeling the dust distribution of the torus. 
In particular, small-scale position-angle-dependent variations of 
the visibility for a given baseline and wavelength seem to be 
clear evidence for the clumpy structure of the torus ( jHonig et al.l 
[2OO6; Schartmann et al. 2008). However, interferometry con- 
strains quite direct information about the radial intensity dis- 
tribution of the torus emission: When observing almost face- 
on-projected type 1 AGN, no significant position-angle change 
of the visibility is expected. Thus, if data for different base- 
line lengths (corresponding to different resolution scales) is 
available, we have a direct and model-independent access to 
the surface b rightness distribution Is of the torus (see also 
iKishimoto et al. 2009a). As shown in Sect. 12.21 Is oc SyOS pri- 
marily depends on a. By modeling the baseline- and wavelength- 
dependence of the visibility, we will be able to tightly constrain 
the radial dust distribution, mostly independent of the SED. 

In Fig. [TT] we show the baseline-dependence of the visibil- 
ity for four different wavelengths. The wavelengths have been 
selected to be in line with current and future interferometric in- 
struments at the VLT and the Keck interferometers: the /T-band 
at 2.2 //m (e.g. VLTI/AMBER), 4.0 //m in the L/M-band region 
(e.g. KI or future VLTI/MATISSE), and 8.5 and 12.5 yum in the 
A^-band (e.g. VLTI/MIDI). For the simulations, we use fixed 
N{) - 7.5 and 7?out = 150 and present results for different ra- 
dial power law indices a - -0.5, -1.0, -1.5, and -2.0. The top 
and bottomrows in Fig.[TT]show ; = 0° and i = 90°, respectively. 
For all wavelengths the variation with position angle in steps of 
10° is shown. Each column in Fig. [TTI represents one particular 
random arrangement of clouds for the given set of torus param- 
eters. 

We will first concentrate on the type 1 case where the ob- 
served visibilities can be directly used to constrain the dust dis- 
tribution. For the shown baseline lengths from 0-450 m, the 
/T-band visibilities in the visibility range ^0.3 have an overall 
similar shape with only little dependence on a. For instance, 
the ^-band visibility at 150 m is 0.5 for a = -2.0 and 0.4 for 
a - -0.5. This similarity can be understood since most of the K- 
band emission is originating from the hot dust which is located in 
the innermost region of the torus, regardless of the actual torus- 
internal obscuration or dust distribution. 



In a smooth dust distribution we would expect no position- 
angle dependence of the visibility in a type 1 torus. The fact 
that we see this PA variations is an evidence for dumpiness. 
The strength of these variations is depending on wavelength. 
Individual dust clouds are optically thick in the /T-band (see 
Sect. 13.1.2b . so that the visibilities can vary between different 
PAs due to the actual clouds' placements and resulting torus- 
internal obscuration effects. At longer wavelengths, the individ- 
ual clouds are optically thin and the variations with position 
angle becomes smaller. The variations still seen at these wave- 
lengths are coming from a combination of "clumping of clumps" 
in the torus, resulting in both stronger local emission and ab- 
sorption, or from "holes" in the distribution. From this analysis, 
we can draw an interesting first conclusion: The strength of the 
position-angle dependence of the visibility in type 1 AGN is re- 
lated to properties of the clouds. This concerns, in particular, the 
optical depth of the clouds which is related to the cloud dust 
mass or geometrical properties and the dust composition. In ad- 
dition, statistical effects will play a role: the same overall torus 
parameters (a, No, Oo), can be kept constant when increasing the 
cloud sizes /?ci by decreasing the total number of clouds in the 
torus (see Sect. 12.21 ). Less number of clouds increase statistical 
effects in the distribution (i.e. holes, obscuration along certain 
line-of-sights) so that the "dumpiness variations" increase. 

As previously mentioned, the /T-band visibility curves at vis- 
ibilities ^0.3 are very similar when comparing different a, so that 
it can be used as a reference for the other wavelengths. As can be 
seen in Fig.[TTl the baseline-dependence of the visibility changes 
with a for longer wavelengths: The A^-band visibilities increase 
from yi2.5/im = 0.15 to 0.65 at 150 m when a changes from -0.5 
to -2.0. The basis of this behavior is the spatial regions and tem- 
peratures involved in emitting light at the given wavelengths. At 
short wavelengths, only hot dust contributes to the emission and 
no light will be emitted from cool dust. This makes the emis- 
sion region small, no matter how the dust is distributed. In a 
simplified picture of a type 1 AGN, we expect that only a very 
small ring is seen, approximately with the size corresponding 
to the equilibrium temperature peaking at this wavelengths. At 
longer wavelengths both hot and cooler clouds contribute to the 
flux. Thus, the emission region appears more extended. The ac- 
tual extension, however, depends on how the dust is radially dis- 
tributed. If the density is steeply decreasing outwards, the size 
is still small and comparable to only hot emission. If the density 
distribution is flat, larger distances contribute significantly and 
the size is much more extended than for hot emission. 

The described principles are also valid for type 2 objects. 
However, in type 2 AGN obscuration within the torus becomes 
more important. In the bottom row of Fig. [TT] we show the 
baseline-dependence of the visibility for an inclination angle / = 
90°, i.e. edge-on as in type 2 AGN. Since clouds at larger radii 
along our line-of-sight obscure clouds at smaller radii (where 
most of the hot emission is originating), the brightness contrast 
between regions close to the AGN and further away decreases, 
so that the tori appear more extended. Thus, the visibilities for 
a given baseline and wavelengths are generally smaller than in 
the type 1 case. Since obscuration plays an important role, the 
small-scale variations of the visibility are stronger. In addition 
and as expected, the tori show elongation along the mid-plane 
which is responsible for a large fraction of the position-angle de- 
pendence of the visibility shown in Fig. [TT] (bottom). Similar to 
the type-1 case, the /T-band visibility shows the strongest small- 
scale visibility variation because the clouds are optically thick. 
The essential dependence of the visibility curves on wavelength 
are also comparable to the type 1 case: Longer wavelengths show 
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Fig. 11. Dependence of the model visibility on the baseline for an AGN at 15Mpc distance and a sublimation radius (- near-IR 
reverberation radius) of 0.025 pc. Each column represents simulation results for a different radial power law index a - -0.5, - 1 .0, 
-1.5, and -2.0, respectively, (top: type 1 AGN at / - 0°; bottom: type 2 AGN at / - 90°) with constant No - 7.5 and Rom - 150. 
In each panel, we show visibility curves at 2.2 (red-dashed lines), 4.0 (green-dotted lines), 8.5 (light-blue solid lines), and 12.5 fim 
(dark-blue solid lines). The different lines per wavelength reflect the position-angle dependence of the visibility in steps of 10° and 
represent one particular random distribution of clouds. 



a stronger change from steep to flat for a changing from -0.5 to 
-2.0. It is interesting to note that at a = -2.0, the type 2 AGN 
visibilities are very similar to type 1 AGN. This can be under- 
stood when considering that in the steep cases, most of the dust 
is located within few sublimation radii. Thus, the overall emis- 
sion is dominated by hot clouds without obscuration from larger 
scales. 

In summary, as illustrated in Fig. [TT] the baseline- 
dependence of the visibility in type 1 AGN directly traces the 
radial dust distribution in the torus: If the dust distribution is 
steep (fl = -1.5 or -2.0), the baseline-dependence of the visi- 
bility is comparably flat, and vice versa, and the effect becomes 
stronger for longer wavelengths. In principle, this also applies 
to type 2 AGN, but obscuration effects and small-scale visibil- 
ity variations make the a-dependence of the visibility less obvi- 
ous. Although we only show results for A^o = 7.5, we note that 
these conclusions are mostly independent of A^o for the range of 
A^o = 2.5 to 10 that we investigate in this paper. 

4.2.2. The wavelength-dependence of the visibility 

In Fig. [TT] we can see that for the "extended" tori with a - 
-0.5 and -1.0, the visibility curves at wavelengths longer than 
2.2 /zm are below the TiT-band visibility curve. In other words, 
the wavelength-dependent size of the torus grows faster than the 
wavelength-dependent resolution decreases for a given baseline. 
On the other hand, the "compact" tori with a = -1.5 and -2.0 
show the opposite behavior, i.e. the long wavelength visibility 
curves are close to or above the /T-band curve for visibilities be- 



low 0.5. Of course, this reflects the fact that most of the dust 
is confined to the inner torus region, making the apparent size 
small at all wavelengths. While this finding is true for both type 
1 and type 2 orientations, it is much more obvious in face-on 
geometries where obscuration effects do not play a significant 
role. 

Instead of looking at the baseline-dependent visibility at 
different wavelengths, we can also compare the wavelength- 
dependence of the visibility for given baselines. In Fig. fT2l we 
show visibilities in the 2 to 40 //m range for both type 1 (top 
row) and type 2 (bottom row) orientations of the torus. As in 
the previous figure, we plot curves for a - -0.5, -1.0, -1.5, 
and -2.0 (from left to right column), A^o = 7.5, /?out = 150, and 
6o =45° are fixed. Each panel contains simulations for baseline 
lengths of 45 m (red), 60 m (green), 90 m (light blue), and 150 m 
(dark blue). For each configuration, various lines illustrate the 
position-angle dependence. The visibility usually drops strongly 
with increasing wavelength up to the mid-IR in the extended 
cases (a = -0.5 and - 1 .0) for both type 1 s and type 2s. The turn- 
around to increasing visibility is depending on baseline and oc- 
curs at ~ 14 yum for 45 m, ~ 18 jum at 60 m, and ~ 27 jum at 90 m. 
This trend is less obvious in the type 2 cases where the 10 and 
18 fim silicate features are seen in the visibilities. The overall in- 
crease in turn-around wavelength for increasing baseline length 
goes along with the increase in spatial resolution. From this we 
conclude that up to the turn-around wavelength, the size of the 
corresponding emission region increases faster than the loss of 
wavelength-dependent spatial resolution for a given baseline. On 
the other hand longer wavelengths would correspond to cooler 
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Fig. 12. Wavelength-dependence of the model visibilities for different baselines for an AGN at 15 Mpc distance and a sublimation 
radius (- near-lR reverberation radius) of 0.025 pc. Each row represents simulation results of one random arrangement of clouds for 
different radial power law index a (top: type 1 AGN at / - 0°; bottom: type 2 AGN at / = 90°) with fixed Nq = 7.5 and /?out = 150. 
In each panel, we show visibility curves for baseline lengths 45 m (red lines), 60 m (green lines), 90 m (light-blue lines) and 150 m 
(dark-blue Unes). The difi^erent lines per wavelength reflect the position-angle dependence of the visibility. 



dust in the outer regions of the torus where both obscuration 
effects suppress direct AGN heating and our model space ends 
(^oiit = 150 corresponds to cloud temperatures of 70 - 180K 
or > 16 yum). Thus, the size of the emission region cannot grow 
any further. In the case of compact dust distributions (a = - 1 .5 
and -2.0) all clouds are confined to small radii so that there are 
almost no clouds contributing from cooler regions. As a result 
the size of the torus is very similar at all wavelengths, so that the 
loss of resolution with increasing wavelength for a fixed baseline 
is dominating the trends of visibilities. 

The difference of the wavelength-dependent visibility for 
different selections of a opens an interesting possibility to nar- 
rowly constrain the radial distribution of the dust, at least in type 

1 AGN: By obtaining interferometric visibilities for an AGN at 

2 different wavelengths for the same baseline length - preferably 
one visibility point in the near-IR, one in the mid-IR - it is pos- 
sible to constrain the radial dust distribution index a or, at least, 
characterize the radial dust distribution as "extended" (a > - 1 .0) 
or "compact" (a < -1 .5). Such experiments can, for example, be 
done using the UT telescope baseline configurations at the VLTI 
in combination with the AMBER and the MIDI instruments. On 
the other hand, since MIDI akeady provides some wavelength 
coverage from 8 to 1 3 jjm, using MIDI alone may already give 
some hint. In Fig.[T3]we show model visibilities for A^o - 2.5 
to 10 for a - 0.0 to -2.0 for a baseline of 60 m. 0o - 45° and 
^out = 150 have been fixed. The strongest difference can be seen 
between the extreme cases in type Is: While overall trend in vis- 
ibilities for a = -2.0 is flat for all A^o (same visibilities at 8 
and 13 i-im), the shallow distributions (e.g. a - -0.5) show some 
clear trend of decreasing visibility in the same wavelength range. 



As shown in Fig.[T3]the silicate feature appears in the visibil- 
ities. But unlike in the SEDs, it is very shallow for small A^o and 
seems in "absorption" for steep dust distributions and in "emis- 
sion" for flat ones in case of type Is. In type 2s it is always in "ab- 
sorption". This tells us that the feature-emitting region is more 
extended with respect to the continuum in steep distributions and 
type 2s while it is more compact in shallow distributions. Given 
the dependence on A^o (larger A^o increases both effects), this 
difference can be ascribed to torus-internal obscuration effects. 
Since optical depth effects along the line-of-sight are responsi- 
ble for the exact strength of the feature, we expect that, aside 
from a and A^o, other parameters will also have an influence, such 
as 00 (as shown in Sect. l4.1.3l for the SED), the dust chemistry 
and the optical depth of the individual clouds. Interestingly, at 
around a - -I the type 1 visibilities do not show the feature 
which might indicate canceling out of some of the parameter de- 
pendencies. The presence of a silicate "absorption" feature in 
the visibilities of type 2 geometries have been claimed to be a 
unique feature to clumpy torus models. Smooth dust distribu- 
tions are allegedly not showing this feature in the visibility spec- 
trum (iJaffe et al.1 120041: ISchartmann et al.1 12005). On the other 
hand, a silicate absorption featur e is promin ently displ ayed in in- 
terferometric data of NGC 1068 dJaffe et al. 2004; Po ncelet et al.1 
l2006l:lRaban et al.l2009}) . Here, we show that clumpy models can 
show a variety of feature appearances in the visibilities for both 
type 1 and type 2 AGN. 
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Fig. 13. Wavelength-dependence of the model visibilities in the MIDI wavelength range from 8-13 jum for an AGN at 15 Mpc dis- 
tance and a sublimation radius (= near-IR reverberation radius) of 0.025 pc at a baseline of 60 m. Each column represents simulation 
results of one random arrangement of clouds for different A^o ranging from 2.5 (left) to 10 (right). In each panel we show models for 
radial power law indices a - 0.0 (red), -0.5 (orange), -1.0 (green), -1.5 (light blue), and -2.0 (dark blue). The top row represents 
a type 1 AGN at / - 0°, the bottom row a 2 AGN at ; = 90°. The different lines per a model reflect the position-angle dependence 
of the visibility. 



4.2.3. The constraining power of IR interferometry on torus 
properties 

As a summary of the model results for IR interferometry we 
emphasize the constraining power that IR interferometry has 
when combined with proper modeling. Here the radial dust dis- 
tribution is particularly accessi ble in type 1 A GN. This has 
been demonstrated recently by iKishimoto et al.l (l2009b.) on a 
set of interferometric data of type 1 AGN and in our pa- 
per I for NGC 3783 ( Honig et all boiOl) . The main reason for 
this constraining power is that IR interferometry can measure 
wavelength-dependent sizes of the continuum emission region 
which corresponds to the temperature distribution inside the 
torus and is, in turn, determined by the dust distribution and the 
source function of the clouds (see eq. |6]l. Since we showed in 
Sect. 14. 1.51 that the mid-IR continuum slope is also determined 
by the radial dust distribution, a correlation between mid-IR 
spectral index a and the compactness in interferometric obser- 
vations can be expected (compact sources have bluer a). 

The silicate feature observed in visibilities has less constrain- 
ing power on the models. The reason is that there exists some de- 
generacy in parameters which control the strength of the feature. 
We have seen in Fig. [13] that the "absorption" and "emission" 
feature in type 1 AGN change with a contrary to what has been 
seen in the SEDs (see Fig.|7]i, which essentially tells us how the 
size of the emission regions of the feature change. The actual 
strength, however, is still influenced by A^o (see Fig. [T3] l and 6q 
(as in Fig.|9]l, and probably also by the optical depth and chem- 
ical composition of the individual dust clouds. It may be still 



possible that simulatneous modeling of SED and interferometry 
can break some of the degeneracy and provide information about 
the dust cloud properties (optical depth) using the strength of the 
"dumpiness variations". 



5. Summary and conclusion 

We presented modeling results from our upgraded radiative 
transfer model of 3D clumpy tori and described how observa- 
tions can constrain some of the parameters used in the model. 
The upgrade concerns the handling of the diffuse raditation field 
within the torus for which we presented a statistical approxi- 
mation. We outlined our modeling strategy and introduced the 
mathematical foundation. The model parameters have been con- 
nected to properties of the torus, such as the volume filling factor 
or the cloud number density. We also derived a simplified type 1 
torus model from the full and more complex model. This simpli- 
fied model has been used to illustrate which geometric parame- 
ters have the strongest impact on the torus emission, which are 
the radial dust cloud distribution power-law index a, the mean 
number of clouds along a radial equatorial line-of-sight A^o, and 
the torus height described either as the scale height h or the half- 
opening angle 6q. We also noted that the cloud radii R^x do not 
have a significant impact on the torus simulation results, while 
the cloud dust properties (ry and dust composition) will proba- 
bly be more important. 

In the next step, we presented results of Monte Carlo 
simulations of dust clouds using different dust compos itions 



We used standard ISM compositions 



;rent gust compos itions, 
with iDraind (l2003i) and 
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lOssenkopf et al.l (1 19921) silicates, respectively, and compared the 
cloud SEDs to simulations with a large grain composition and 
Graphite-dominated dust. The differences between the cloud 
SEDs from different dust compositions are an effect of the opti- 
cal depth in the infrared and the different (Planck mean) opaci- 
ties of each dust distribution which result in a slightly different 
temperature gradient with distance from the AGN. 

By using the pre-calculated cloud SEDs and images, we sim- 
ulated torus SEDs and studied how different parameters affect 
the model SEDs. The SEDs are presented in a suitable way to be 
scaled to observations. Our main conclusions are: 

• The dust clouds with standard ISM dust seem to provide 
reasonable torus SEDs which seemingly capture the essence 
of recently observed IR characteristics of AGN (e.g. moder- 
ate lOyum silicate emission and absorption features). While 
not showing explicit dependence of the SEDs on the optical 
depth of the individual clouds, t^i = 30 - 80 for ISM stan- 
dard dust is probably the range that provides good agreement 
with observations. 

• The torus SEDs, in particular in type 1 cases, are quite sen- 
sitive to the selection of a. Distributing more dust clouds at 
larger distances (e.g. a - -0.5 or -1.0) results in redder 
mid-lR colors. We showed that a is the main parameters that 
determines the mid-lR spectral index a with only little de- 
pendence on A^o or Oo- Therefore it is possible to use suitable 
mid-lR data to directly constrain the radial distribution of 
the dust in the torus. A formal correlation function between 
a and a has been extracted from the models (see Sects. l4m 
& I4T5] and Fig. [Toll. 

• The appearance of the silicate feature - either in absorp- 
tion or emission - is not only viewing-angle (inclination) 
dependent, as could be naively concluded, but depends on 
the radial and vertical distribution of the dust (i.e. a and Oq 
or h), the average number of clouds along an equatorial line 
of sight (A^o) and possibly on and the dust composition 
(see Sects. 14.1.21 & I4.1.3I I. In fact, all the cloud and torus 
parameters influence the exact shape and depth of the sili- 
cate feature towards more absorption or emission and it may 
be difficult constraining individual parameters only based on 
silicate feature observations. That being said it is, however, 
crucial that parameters fitting the overall SEDs of observa- 
tions are also in agreement with the silicate feature strength. 
We showed that compact dust distributions (a = -1.5 and 
-2.0) with Tci = 50 do not produce any silicate absorption 
features for A^o ^ 10 in type 2 line-of-sights. Moreover, the 
silicate emission features in type 1 AGN become more pro- 
nounced with decreasing A^o and/or more compact a. From 
the commonly seen weak silicate emission features in type 
1 AGN and silicate absorption features in type 2 AGN, we 
would expect that, in general, AGN dust tori should have, on 
average, rather extended dust distributions with No ^ 5. 

• We present an easy way to compare torus models to observa- 
tions. For that we plot the apparent silicate feature strength 
Tsi against the mid-IR spectral index a which are the main 
observational properties of torus spectra (see Sect. l4.1.3] l. On 
top of this, we overplot the parameter range that is covered 
by the models. Neglecting inclination effects, a is almost ex- 
clusively depending on a while Tsi is influenced by A^o, ^^o, 
and the cloud dust properties (ry and the dust composition). 
If we select a separate sample of type 1 and type 2 AGN 
(i.e. having average inclinations of / ~ 30° and / ~ 75° re- 
spectively) and consider sample average type-l/type-2 ratios 
(leading to a an average do ~ 40 - 50°) it should be pos- 



sible to provide some sample average constraint on a com- 
bined No-Tci ■ Considering both 1 //m and 1 8 jum silicate fea- 
ture s might help breaki ng some of the remaining degeneracy 
dThompson et al.ll2009h . The principle of this approach ha s 
been introduced in paper I of this series jHonig et al.ll2010l) . 

• When modeling near- to mid-IR data, the outer radius does 
not have a significant influence on the near- and mid-IR 
emission for steeper dust distribution than a ^ -1. In these 
cases the majority of the dust is confined to distances <s: Rom. 
On the other hand, for a <: - 1 the dust distribution function 
does not converge for /?out — > so that /?out influences the 
peak of the overall torus emission. It is therefore important 
to properly select 7?out in all cases to avoid artificial cutoffs in 
the brightness distribution of the targeted wavelength range 
(see Sect l4.1.4l i. 

From torus model images, we calculated visibilities to simu- 
late observations with interferometers. In addition to the SEDs, 
these data provide spatial information at different wavelengths 
and at different scales depending on the baseline length. We con- 
clude: 

• In type 1 AGN, the baseline-dependence of the visibility 
directly traces the radial distribution of the dust (here: the 
power law index a) without significant influence from other 
parameters (see Sect. 14.2.1b . The general principle is also 
valid for type 2 AGN, but obscuration and geometric effects 
result in significant dependence of the visibility on position 
angle. 

• By comparing baseline-dependent visibilities at two different 
wavelengths (preferably one in the near-IR, one in the mid- 
IR) for a given baseline, one can directly conclude if the dust 
distribution is compact (a > -1.5) or extended (a < -1.0). 
This can be done for both type 1 and type 2 AGN, although 
type 2s are less constraining due to obscuration effects (see 
Sects. |4TT]&|432]|. 

• Small-scale position-angle variations of the visibility in type 

1 AGN are usually associated with dumpiness of the torus 
(see Sect. 14.2.1b . The strength of these variations is con- 
trolled by (1) the optical depth of the cloud and, thus, 
also the dust composition (they are strongest if » 1 at 
the observed wavelength A), (2) the cloud size R^i, and (3) 
the number of clouds in the torus (i.e. statistics which are re- 
lated to the choice of A^o)- It may be difficult to constrain any 
of these parameters individually by observing the variations. 

• The visibilities show silicate features in both type 1 and type 

2 orientations. For type 1 AGN the appearance in "absorp- 
tion" or "emission" has some strong dependence on a and 

the feature is mostly canceled out around a 1 . However, 

the actual strength of the feature is degenerate with several 
parameters as in the case of SEDs (see Sects. l4~Z2l & l4.2.3l l. 
Since the silicate feature in the visibilities is the result of 
both changes of flux and emission region size, it cannot be 
directly associated with emission properties of or within a 
certain region. 

The presented models can be used as direct input for 
interpreting observational data. We are pleased to provide 
model SEDs and images based on our clumpy torus model 
{"CAT3D") to interested colleagues on request and via the web- 
page |http : //cat 3d . sungrazer . org 
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Appendix A: AGN radiation 

As AGN/accretion disk input spectrum, we use a four- 
component power-law which is approximately in agreement 
with the range of power-law fits to QS O composite spectra (e.g . 



Zheng et al 



Scott etalJlioCT : 



1997UManske et al.lll998tlVanden Berk et"ani20QU 



.1/3 



..-2 



V < 10'"^ Hz <-> > 3;um 
lOi"*- 10'^ Hz ^ 0.3-3;um 
lO'^-lO'^'Hz ^ 0.03 -0.3 yum 
> 10"= Hz <^ < 0.03 //m (> 0.4 keV) 



(A.l) 



The main difference to the spectrum used in lHonig et al.l (l2006l) 
concerns the relatively blue part with Fy oc v'''^ in the op- 
tical to near-IR between 0.3 and 3 fim. This spectral shape 
has recently been confirmed by near-IR polarimetry of quasars 
(lKishimotoetal.l l2008 ) . 

We note that the actual shape of the spectrum is much less 
important for the heating of the dust than the total (spectrally- 
integrated) luminosity of the accretion disk, Lboi - J Lydv. 
In the optically-thick case, the total absorbed heating power is 
^abs - j Qabs(v)Lv(iv ~ 0.5 X Lboh depending on the actual dust 
composition, where 2abs(v) is the absorption efficiency of the 
dust. Thus, the spectral shape of the incident radiation is inte- 
grated out in the process of radiative transfer. 



Appendix B: The influence of the accretion disk on 
IR interferometry of dust tori 

The presented torus models do not include any accretion disk 
contribution. For type 2 AGN, it is safe to conclude that most of 
the observed IR flux is coming from the dust torus since, by def- 
inition, the accretion disk has to be obscured. On the other hand, 
depending on wavelength, the accretion disk might contribute to 
flie overall AGN IR flux in type 1 AGN. 

Recently, we presented an analysis of the near-IR emission 
of nearby type 1 AGN (Kishimoto et al. 2007). We used a color- 
color diagram to disentangle the hot dust emission from the ac- 
cretion disk power law spectrum extending from optical wave- 
lengths into the IR. A simultaneous power-law+black-body fit 
suggested that the accretion disk contribution to the AGN near- 
IR is around 10-50% in the //-band and only 5-25% in the 
A'-band. The decrease of accretion disk contribution with wave- 
length is consistent with the recent obs ervation of the blue tail 
of the near-IR accretion disk emission dKishimoto et alJl2008l) 
which suggests that for any wavelength longer than 2 yum, the 
ac cretion disk has margin al effect on SED and interferometry. 
In iKishimoto et al] (l2007h . we also showed that a possible 20% 
accretion disk contribution in the /T-band increases the visibil- 
ities by at most $0. 1 for the longest baseline shown in Fig [121 
Since the position-angle variations in the near-IR are of the same 
magnitude, any accretion disk contribution will be hidden under 
these variations. Thus, we refrained from adding the accretion 
disk contribution to the presented figures and note that for inter- 
pretation of near-IR interferometry, a case-to-case analysis based 
on the color-color diagram should constrain to what degree the 
accretion disk has to be taken into account. 
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